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Abstract. Using the Anderson model for disordered systems the fluctua- 
tions in electron spectra near the metal-insulator transition were numer- 
ically calculated for lattices of sizes up to 28 x 28 x 28 sites. The results 
show a finite-size scaling of both the level spacing distribution and the vari- 
ance of number of states in a given energy interval, that allows to locate 
the critical point and to determine the critical exponent of the localization 
length. 



The statistical description of energy spectra of disordered quantum sys- 
tems is based on the random-matrix theory [1-3]. One of the main prop- 
erties of the spectra of random-matrices is the repulsion between their 
eigenvalues. For a disordered metal, such a correlation of energy levels is 
caused by a pronounced overlap of delocalized one-electron states. By in- 
creasing the disorder of a random potential the system is known to undergo 
a metal-insulator transition (MIT). On the insulating side of the MIT, the 
energy levels of the localized states are not correlated due to vanishing of 
the level repulsion. A central question in the problem of level statistics is 
how a character of spectral fluctuations varies when the system changes 
from the delocalized to the localized regime. 

On the metallic side of the MIT the distribution function P(s) of neigh- 
boring spacings between levels can be described by the Wigner formula [4] 

p w{s) = - s exp ( - -s 2 \ , (1) 

where s is measured in units of the mean inter-level spacing A. In the 
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localized regime the spacings are distributed according to a Poisson law 

P P (s) = exp(-s). (2) 

The crossover of P(s) between the Wigner (1) and the Poisson statistics (2) 
which accompanies the MIT was extensively investigated both numerically 
[5-8] and analytically [9, 10]. It was pointed out by Shklovskii et al [7] 
that the level spacing distribution exhibits critical behavior near the MIT. 
Therefore we focus on the finite-size scaling properties of the distribution 
P(s). In addition, we analyze the variance of the number of energy levels 
([5N(E) 2 ]) as a function of the average level number (N(E)) in a specified 
energy interval E, and study its universal peculiarities at the critical point. 

One of the simplest models describing a disorder-induced MIT is the 
Anderson model 

H = £ n a n «n + ( a n a m + ^4)- (3) 

Here a+ and a n are the creation and annihilation operators of an electron 
at a site n in a lattice; m denote the nearest neighbors of n. The on-site 
energy e n is measured in units of the overlap integral between adjacent 
sites and is uniformly distributed variable in the interval from —W/2 to 
W/2. The parameter W specifies the degree of the disorder of the system. 
The critical disorder of the MIT which occurs in the middle of the band 
corresponds to W c ~ 16.5 [11]. For a finite system when W <C W c , the level 
statistic is close to (1) and when W S> W c it obeys (2) [5]. 

In order to find the electron spectrum in the critical region we diag- 
onalized numerically the real symmetric Hamiltonian (3) for simple cubic 
lattices of the size LxLxL with periodic boundary conditions. The Lanczos 
algorithm for eigenvalue problem was used in a version that was especially 
designed for very sparse and big matrices with hierarchic structure. We 
applied the algorithm for cubes with L = 6, 8, 12, 16, 20, 24, 28 at various 
degree of the disorder close to W c . We consider an energy interval which is 
centered at e = and has such a width that it contains a half of all eigenval- 
ues. As we are interested in sample-to-sample fluctuations in the spectra 
the calculations were carried out for ensembles of different random config- 
urations. After unfolding the spectrum the histograms of several spectral 
distributions were constructed by use of 10 5 spacings calculated for each 
pair of {L, W}. 

Fig. 1 displays the distribution function P(s) calculated near the MIT 
for two cubes of sizes L = 6 and L = 28. By increasing the disorder W the 
spacing distribution for both L varies continuously from P\y{s) to Pp(s) 
over all range of spacings. Results of calculations for other, intermediate 
sizes confirm this continuous crossover. But the behavior of this crossover 
substantially depends on the size of the cube. One can see that P(s) changes 
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Figure 1. The level-spacing distribution P(s) for various disorders near the transition. 
Continuous curves correspond the Wigner (1) and the Poisson (2) distributions for the 
metallic and insulating phases, respectively. 



faster between (1) and (2) for L = 28 than for L = 6. The size-dependence 
of P(s) is observed on both sides of the MIT. However at the transition 
point when W c = 16.5 the spacing distribution has almost the same form 
for all L from 6 to 28. The independence of the level spacing distribution on 
the size of system at W c is in good agreement with predictions of the new 
universal level statistics at L — > oo which exists exactly at the transition [7]. 

In order to study the scaling properties of the spacing distribution in 
more detail we introduce the quantity a = Jq(P(s) — Pw(s))ds/ Jq(Pp(s) — 
Pw(s))ds, which describes a normalized deviation of P(s) from the Wigner 
distribution (1). We chose s = 0.473, the crossing point of Pw(s) and 
Pp(s), in order to study P(s) in the range of small spacings. In this case 
a corresponds to a "strength" of the repulsion of two consecutive levels 
when the separation between them is less than average spacing A. In the 
thermodynamic limit a = for W < W c , and a = 1 for W > W c . For finite 
L it is reasonable to assume a scaling law 

a(W,L) = f(L/£(W)), (4) 

where £(W) is the correlation length of the transition. 

Fig. 2 shows the L-dependence of the parameter a near the critical 
point. In the insulating regime a grows with L approaching its limit a = 1, 
whereas in the metallic phase a decreases to zero. The change of the sign of 
the size effect takes place at the critical point. Therefore we can determine 
very accurately the critical value of the disorder W c = 16.35 ±0.15. Similar 
results for smaller systems were earlier obtained in Ref. [7, 8] for so = 
2 which corresponds to using P(s) only in the range of asymptotically 
large spacings. In order to obtain the scaling curve (4), on which all data 
points collapse, we shifted InL by ln£(W) separately for each W. The 
overlap between adjacent values of W allows to fit most of the points onto 
two branches for W < W c and W > W c corresponding to localized and 
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Figure 2. The quantity a as a function of Figure 3. One-parameter scaling be- 
the size L at different degrees of disorder near havior of the level statistics. The quan- 
the metal-insulator transition. tity a as a function of L/£(W). 



delocalized regimes, respectively (Fig. 3). As s result one finds the disorder 
dependence of the correlation length i(W). We calculated also the critical 
exponent of the localization length v. Using the singularity of £ near the 
transition point (,(W) oc| W — W c \~ v one can expand the relation (4) 
to a power series. Taking into account only the linear term, a(W, L) = 
a c (L) + const L l / U (W — W c ), and applying the x 2_ criterion to fit the data 
plotted in Fig. 2 we found that v = 1.45 ± 0.1. 

Another important quantity which is used to describe the level statistics 
is the magnitude of fluctuations of the number of energy levels 5N(E) 
in a given energy interval E. The variance ([5N(E)] 2 ) characterizes the 
"stiffness" of the spectrum. From the statistical viewpoint it is reasonable 
to investigate the dependence of the variance ([5N(E)] 2 } on the average 
number of levels (N(E)) in the vicinity of the MIT. Here (...) denotes the 
averaging over the random configurations with the same disorder W. In the 
metallic regime the variance is known to be defined by the Dyson formula, 
(5 2 N) M = 2/tt 2 (ki{N)+C), where C » 2.18, provided that (TV) » 1 [1]. In 
the strongly localized regime the levels are not correlated, hence {5 2 N)\ = 
(N), that is much larger than {8 2 N)m- Decreasing disorder suppresses the 
fluctuations SN, so that the variance changes from {[5N] 2 )i to {[SN} 2 ) M [5]. 
Exactly at the transition the variance behaves linearly with the average level 
number, (5 2 N) C = k(N) [5, 10], as in the insulating regime. However the 
numerical factor k is less than unity. 

We calculated the dependence of the ratio ([5N(E)] 2 ) / (N(E)) on the 
average number of levels within a given interval E for different lattice sizes 
L at the disorder W varying from 12 to 20, as shown in Fig. 4. It is obvious 
that near the MIT this ratio exhibits the critical behavior. In approaching 
the transition, W = W c , it becomes size-invariant (Fig. 4c). Such a behavior 
is closely related to the universality of the critical level spacing distribution, 
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Figure 4- The variance {[5N(E)] 2 } as a function of (N) for various sizes L at disorder 
W: a) 12 and 20; b) 15 and 18; c) 16 and 17 for lower and upper set, respectively. Critical 
variance at W c = 16.5 is shown for different L: a) 24 (+); b) 28 (x). Solid line corresponds 
to the Dyson law. 



as mentioned above. In addition, one can observe that at the MIT the 
relative fluctuations of the level number ([5N(E)] 2 ) / (N(E)) decreases with 
the energy, when (N) = E/A ~ 1, and then varies very weakly over two 
orders of magnitude of (N), tending to a constant value k ss 0.32. For 
E > 200A the numerical errors becomes larger due to the finite number 
of realizations. Our results are consistent with the suggestion about the 
proportionality between the variance and the average level number [5], but 
deviate from the power law proposed recently in [13]. 

Determining the ratio (5N(E) 2 )/(N(E)) as the function F(E,L,W), 
one can analyze the scaling properties similarly to those of P(s). For ex- 
ample, the disorder dependence for the given energy width E = 20A is 
shown in Fig. 5. Near the critical point, \W — W c \ < 1, the function can 
be linearized F(E,L,W) = k + AL 1 ^ (W - W c ). The factor A depends 
on the energy E, in contrast to that for a, whose critical behavior is not 
sensitive to the choice of s. In order to obtain the one-parameter scaling 
law one should take in account the typical volume per one electron state 
lying in the interval E, L 3 = L 3 (A/E), instead of the total volume L 3 , i.e. 
F -hoc (L /i) l / u . Fig. 6 shows {5N 2 )/{N) vs. L /£(W), where £(W) is 
taken from the analysis of a. All data belong to the common two-branch 
curve regardless to L, W, and chosen {N). Indeed, the fluctuations 5N in 
the interval E are mainly defined by the states confined within the correla- 
tion volume £, which are separated by the energy = A (L/£) 3 . The num- 
ber of those states is N 6 = E/A 6 . Therefore F = k + N^ 1/3u sign(W -W c ). 

In conclusion, we have used the level spacing distribution P(s) as a scal- 
ing variable in order to detect the critical behavior at the disorder-induced 
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Figure 5. Disorder dependence of the ra- 
tio (SN 2 )/{N) at the given number of lev- 
els (N) = E/A = 20 for various sizes L. 




Figure 6. {5N 2 )/(N) function of 

Lo/£{W) for different L and various (N) 
from 10 to 100. 



metal-insulator transition. The results for the critical exponent v and the 
critical disorder W c , which were obtained by numerically diagonalizing the 
Anderson Hamiltonian for up to 28 3 lattice sizes, are consistent with those 
obtained earlier by using completely different approaches [12]. Our calcu- 
lations showed that the universality of the level statistics at the transition 
is revealed not only in the form of P(s) but also in the dependence of the 
variance of the number of electron states in an energy interval of a given 
width on their mean number (N). We also analyzed the scaling behavior 
of the function (5 2 N) = F((N)). At the MIT this function is found to be 
L-invariant with a leading linear term, (5 2 N) C = k(N), where k « 0.3. 
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